Can infrared gravitons screen A? 
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It has been suggested that infrared gravitons in de Sitter space may lead to a 
secular screening of the effective cosmological constant. This seems to clash with 
the naive expectation that the curvature scalar should stay constant due to the 
Heisenberg equation of motion. Here, we show that the tadpole correction to the 
\ local expansion rate, which has been used in earlier analyses as an indicator of a 

' decaying effective A, is not gauge invariant. On the other hand, we construct a 

(N ■ gauge invariant operator which measures the renormalized curvature scalar smeared 

over an arbitrary window function, and we find that there is no secular screening of 
. this quantity (to any given order in perturbation theory). 

I. INTRODUCTION 

p • In de Sitter space, long wavelength gravitons are frozen in, producing a cumulative de- 

O \ formation of spacetime on large scales. This can be seen, for instance, in the tree level 
' graviton two point function. In the transverse traceless gauge, the behaviour of gravitons 



is similar to that of massless minimally coupled scalars pLl], and their two point function 
^ ■ {h{x)h{x')) grows logarithmically with scale. Globally, such increasing departure from a de 

•jH ! Sitter metric cannot be undone by a gauge transformation. Nevertheless, infrared gravitons 
do not contribute to tidal forces on small scales. The tree level two point function for the 
O ' Riemann tensor {R{x)R{x')) is infrared finite, and the contribution of gravitons with wave- 
\Q • length much larger than the curvature scale is in fact negligible. Hence, to lowest order 
O I in perturbation theory, the local geometry remains everywhere close to the unperturbed de 
Sitter space ^. 

It has long been suggested that graviton interactions may dramatically alter this picture, 
■ potentially leading to infrared screening of the cosmological constant [2|. The basic idea is 
rS I the following. Gravitons carry energy and hence they are a source of the gravitational field. 
I Hence, it is conceivable that the accumulation of infrared modes crossing the horizon in the 
expanding de Sitter phase may backreact on the average expansion rate of the universe. A 
priori, it is unclear whether infrared gravitons can have much of an effect, since the "energy" 
in the gravitational field is contained in derivatives of the metric. To make a quantitative 
estimate, the authors of jl] calculated the graviton tadpole {hfj_^) at the two loop order. 
From that, they obtained the "tadpole corrected" expansion rate of the universe H{{h^J)), 
which turned out to decrease quadratically with cosmic time. This slowing down of the 
expansion rate of the "averaged metric" , was interpreted as a secular screening of A by the 



(N 



o 



^ Provided, of course that H is well below the Planck scale. The graviton power spectrum is scale invariant 
for wavelengths above H^^ , with amplitude of order h ^ H/Mp. 



long wavelength modes. If true, this would be a spectacular effect of low energy quantum 
gravity, with implications for the cosmological constant problem 

The purpose of this paper is to reanalyze this problem, with an emphasis on gauge 
invariance. In Section II we show that the tadpole correction to the expansion rate, as defined 
in Ref. is not gauge invariant (and can in fact be given an arbitrary time dependence). 
In Section III we discuss a physically motivated gauge invariant definition of the expansion 
rate, which in the present context essentially links it to the local value of the Ricci scalar. 
In Section IV we calculate a gauge invariant smeared expectation value of the Ricci scalar, 
suitably renormalized, showing that there is no infrared secular screening of this quantity. 
Our conclusions are summarized in Section IV. 



II. ON THE TADPOLE CORRECTION TO THE EXPANSION RATE 

The theory under consideration is pure gravity with a cosmological constant. The action 
is given by 

Ssr = ^J v^(7^-2A)^/^x, (1) 

where 71 is the Ricci scalar and k = SnG. Here, G is Newton's constant. We are interested 
in perturbations around the de Sitter space solution, and for definiteness we shall adopt the 
flat chart description. The perturbed metric can be written as 

g^^ix) = a'^i'n)['nf^u + • (2) 

Here 1]^,^ is the Minkowski metric, a{ri) = —l/{Hori), with — oo < r/ < the conformal time 
and Hq the constant unperturbed expansion rate. 

To perform a systematic perturbative expansion using the path integral, we must add 
to ([T]) a gauge fixing term Sgf = — {1/2)1]'^'^ Ffj^[h]F^[h], where the function is such that 
Ff^[h] = selects one representative out of a given gauge orbit. In Ref. jl] this function was 
chosen as 

F,[h]^a(^h;^^-h„+2h';^^ . (3) 

Here, indices are raised and lowered with the Minkowski metric r]^,^. Suitable counterterms 
will be also be needed in order to remove divergences.^ The total action takes the form 

Slot = Sgr + Sgf + Spp + Set , (4) 

where Spp indicates the Faddeev-Popov (FP) ghost terms and Set the counterterms. The 
graviton tadpole is defined by 

( V)^ = / ^^"^ K-^ e^^*°*['^^l e-^^'°*[^'] , (5) 

JcTP 

where the subindex F refers to the gauge fixing function and ip indicates the set of 
dynamical variables: metric perturbations hfj^^y and the FP ghosts and anti-ghosts. The 



^ General Relativity is non-renormalizable, and the number of counterterms needed in Set increases with 
the number of loops at which we calculate our observables. However, the number is finite at any order, 
as it is usually the case in effective field theories. 



closed time path (CTP) version of the path integral is indicated, since we are interested in 
expectation values (rather than matrix elements between in and out vacua). 

The left hand side of Eq. can be computed diagrammatically order by order in pertur- 
bation theory. In the gauge ([3]), the propagator is infrared divergent in the limit of infinite 
volume [2I, so it is convenient to compactify the spatial directions and consider a finite (al- 
though in principle arbitrary large) co-moving volume. If we choose a spatially homogeneous 
initial state, symmetry requires that 



(6) 



where, = (dr,)^. 

From ([2]) and ([6]), the "averaged" metric {g^v)F is a fiat Friedmann- Robertson- Walker 
(FRW) metric, with expansion rate given by 



HF = H{{h^,)F) 



Ho 



d\n.[a{l + Apf^] 
ail + Ap- BpY'^dri ~ Ap - Bpfl^ 2 (1 + Ap) _ 



1_L ^^'^ 



(7) 



where the prime indicates derivative with respect to rj. 

In Ref. {2!, Tsamis and Woodard calculated Ap and Bp at the two loop order. Upon 
substitution in ([7]) they obtained 



Hp = Hn 



Ho 
in 



6 



+ 0(^6) 



(8) 



which decreases quadratically with cosmological time t = —H~^\n{Hori). As mentioned in 
the introduction, this result was interpreted in as a secular screening of the effective 
cosmological constant by the infrared gravitons. However, as we shall now discuss. Hp is 
not invariant under generic gauge transformations, and so the above interpretation seems 
rather questionable. 

Let us consider a new gauge fixing function G[K\ in the vicinity of F[h]. \{ F[h] =0, then 
we can find a new metric perturbation 



related to h^y by a gauge transformation 5x^ = x'^ 



X^, such that G[h] = 0. Here 



where is the covariant derivative with respect to the full metric g^^ and Xu = duxX^- 
Note that the gauge transformation will in general depend on hf^^, 



(9) 



and even for simple changes of the gauge function F[h] the dependence of x^ h can be 
quite non-trivial. The point, however, is that for every h this transformation will exist. In 
what follows, we shall consider the class of gauge function G in the neighborhood of F which 
are defined through the equation 



F[h] = G[h + 6^h] , 



for some x- 

Both Sgr and are gauge invariant, Sgr[hb\ = Sgr[h + S^h],and Sct[h] = Sct[h + S^h]. 
Moreover {Sgf)p[h] = {Sgf)c[h + S^h] and {SFp)F[h] = {SFp)G[h + S^h]. It is then straight- 
forward to show, by changing variables in ([5]), that 

{h^,u)G = {hf,u + S^h^^)F . (10) 
The variation of the tadpole under gauge transformations is thus given by 

= {[Vxu + h^uhi + [^V + ^A^lxi + h^,,xx^ - '^[v^, + + 0{x^) . (11) 

Here 6^A and 6^B represent the changes in A{ri) and B{ri) defined in ©. In the above 
expression x is treated as a small quantity, but h^i, is not necessarily small. 

In the case when x'^ is a c-number (and by this we mean a function independent of h), 
the only transformation compatible with the symmetries of a flat FRW is 

= f{v)t' ■ (12) 

If we neglect hf^i, in Eq. ( ITTil . the expansion rate H is invariant under gauge transformations. 
The reason is that for a fiat FRW the expansion rate is given in terms of the temporal 
component of the Einstein tensor 

H' = . (13) 

The background is such that 

g;; - A 5^ = . (14) 



It follows that ^x^j^ ~ 0' from (11311 S^H = 0. More explicitly, from ( ITTI) we have 
= —(2////) and S^B = —2/', and linearizing (JTj) we have ^ 

S^H = Hg-Hf = {Ho/2)[6xB - S^A - riid^A)'] = . (15) 

Provided that is independent of h, the above consideration can be extended to the case 
when the tadpole is non- vanishing. Using ([6]) in (fTTi) we find S^A = — (2//?7)(l + A) + A'f 
and S^B = -2/'(l + A - B) - {2f /ri)B + B'f. Substituting these variations in (CD, it is 
straightforward to check that 

Wv) = ^m + --- . (16) 

Here, H{r]) represents the right hand side of (JTj), which depends on time through A and B, 
and the ellipsis denote higher orders in x- 

The simple form of f|T6|) is easily understood. To lowest order in x, the variation S^h^^ 
is exactly linear in /i^,^ [the metric perturbation h^,^ is not treated as a small parameter 



^ A derivation of (|15p under similar assumptions was given in p]. This, however, does not estabhsh that 
H{{hfj,,y)) will be invariant. As shown below, Eq. (flS]) does not hold for generic gauge transformations. 



in Eq. (|TT|) ]. Because of that, {h^J) transforms like a classical metric under infinitesimal 
c-number gauge transformations x^i 

(5^V) = '^x(V) + -- - ' (17) 

Eq. (I16p follows immediately by noting that = (1/3) Gq is a scalar under redefinitions 
of the time coordinate (because it has mixed temporal indices). The gauge dependence (fTBj) 
is therefore rather irrelevant: it indicates that we have changed the parametrization of a 
time dependent function, but this does not change the value of the expansion rate if as a 
function of proper time t [defined by dt = a(l + A — BY^'^dr]]. 

However, the above conclusions do not apply to the case where the gauge transformations 
depend on the metric [1] ^ 

X^ = Xnh]- (20) 
For generic choices of x[h], we should expect 

and from Eq. ffTTj) . we should likewise expect that 

((^xV) ^ '^(x>(V) + ■•• • 

We may thus anticipate that, in general, the expectation value of the gauge transformed 
metric perturbation {h + 6^h) will not be gauge equivalent to the original one (h). 

The expansion rate H{{h^^)) as a function of proper time t would be gauge invariant (and 
therefore meaningful) if and only if {h + 6^h) is related to (h) by a time reparametrization 
(see e.g. the discussion around Eq. ffTTj) ). In equations, this means that for each x'^[h] we 
should be able to find a vector ^'^ such that 

(5xV) = '^?(V) + --- . (21) 

where, by symmetry, C,^ = g{r])t'^. However, an equation of this sort cannot hold for a generic 
x[h]. To illustrate the point, we may restrict ourselves to lowest order in perturbation theory. 



For illustration, note that even the simplest changes in the gauge function F may lead to a complicated 
dependence of x on /i ([^D]) . Consider for instance the one parameter class of gauges 

out of which ([3]) corresponds to a = —1/2. A change of a corresponds to 6F^ — a h,^ Sa, and we have 
By introducing the explicit expression of F and S^h^^,, this equation takes the form 



O^Ah]x'' = h,^ , (19) 

where 0^v[h] is a second order differential operator whose coefficients depend on /i^^ and its first and 
second derivatives. Eq. p9p should in principle be solved in order to find x iii terms of h. It is clear that 
in this case the dependence will be highly non-local (and difficult to find explicitly) , but the point is that 
we cannot restrict consideration to c-number gauge transformations, since generically x depends on h. 



where the expectation value of odd functions of h will vanish. Let us therefore take x to 
be an odd function of the metric perturbation h. In this case Eq. f l2ip to lowest order in 
perturbation theory reads 

{2hxi^X% + VaX^ - ^ = - '^9't^.U • (22) 

Here we have used = g{ri)t^, as dictated by symmetry. When we consider rescaling given 
by x^[h] —>■ k{ri)x'^[h] with an arbitrary function k, we have the terms with k'{7]) and those 
with k{r]) on the left hand side of Eq. ([22]) • In order to satisfy this equality for an arbitrary 
function k, the right hand side also has to have the terms with k'{ri) and those with k{ri). 
This requires that g{ri) should also transform to k{T])g{'r]). Comparing the term proportional 
to k'{r]), we immediately find 

g = {hoxx') ■ (23) 

Then, if we choose = x'^i^At'^, where x° is an arbitrary odd function of /i^j/ (including 
possible arbitrary explicit dependence on rj and contraction of internal indices with the 
vector t^), the = z/ = component of ([221) reads 

(^00,0 X°) = . (24) 

Clearly, this equation does not hold in generic gauges and for generic choices of x^- In 
particular, it does not hold in the gauge defined in Eq. ([3]), which completes our proof ^. 

It follows that the "tadpole corrected" expansion rate is not physically meaningful. 
Rather, given the enormous freedom in choosing x'^ (which can include arbitrary functions 
of h and rj) it appears that H{{h^y)) can in fact be given arbitrary dependence on proper 
time t. 



III. OBSERVABLES 

Gravitational radiation of wavelength shorter than the Hubble radius has an impact 
on the background expansion rate. Formally, this can be accounted for in the so-called 
Isaacson approximation (see e.g. j^), where the Einstein tensor is split into a background 
"long wavelength" contribution and the contribution from short wavelength gravitational 
waves. A classical bath of short wavelength gravity waves does modify the time evolution 
of the scale factor, much like a bath of radiation would. However, this does not mean that 
it screens the cosmological constant, which also contributes to the expansion rate as usual. 

On the other hand, here we are not interested in the effect of short wavelength modes but 
in the collective effect of infrared graviton modes and their interactions. Could these cause 
a secular screening of the cosmological constant? What we mean by this is an adiabatic 
erosion of the expansion rate, such as the one suggested by Eq. ([HI), which would lead to an 
initial quasi-de Sitter phase with 

\H\ < H'^ . (25) 



^ Eq. (|M1) may accidentally hold in some gauge for all possible functions x°- For instance, it holds in the 
transverse and traceless gauge, since (/iqo) vanishes to lowest order. In that case we should examine the 
other components of Eq. (1^^ to see whether they can hold for generic x- 



As emphasized in the previous subsection, instead of calculating the expectation value of 
the metric, it is important to look for some gauge invariant characterization of the expansion 
rate. One such observable was suggested by Abramo and Woodard The basic idea is to 
consider the value of a scalar field (p conformally coupled to gravity, with a constant source 
term J: 

n 



□ + 
6 



= J . (26) 

In a flat FRW universe, the Ricci scalar is given by 

7^ = l2H^{t) + QH{t) . (27) 

For instance, when the scale factor takes the form a oc t^, we have H{t) = pt~^ and TZ = 
(12p^ — 6p)t~'^. The general solution of Eq. (1261) in this case is of the form 

^W = A,«-+A_r-+(^^^jJ|^)^. (28) 

where A± are arbitrary constants and 



—3p ± a/p^ + 4p 

a± = . 

^ 2 

The last term in Eq. (1261) is proportional to t^. As we have ^a± < 2 for p > —4/3, the two 
first terms in fl2^ decay faster than the last term. Therefore at late times only the third 
term will be important. The limit of quasi de Sitter expansion {H ^ H"^) corresponds to 
p ^ 1, and in this case we have 

'^^^^"^2^^) (29) 

We may interpret this result in the following way. The source J creates the field (f). In 
turn, the field is diluted by the expansion, which causes the amplitude to fall off as the 
inverse of the scale factor. At late times, the value of is dominated by the field created 
by the source J during the last expansion time, on the surface of a sphere of radius H~^, 
while the initial conditions set by the coefficients A± become irrelevant. Thus, the late time 
asymptotic behaviour of 0()f:) sourced by a constant J is a measure of the time dependent 
expansion rate H{t), through Eq. f l29l) . Let us now consider the quantum theory. To avoid 
the introduction of a new quantum scalar field in the system, it was suggested in Ref. [3] 
that (pit) be defined as the expectation value of inverse of the retarded conformal propagator 
acting on the constant source 

*W = <^i>. (30) 

This quantity would then be a characterization of the inverse of the square of the expansion 
rate, through the identification (1291) . Although this seems to be an appropriate definition, it 
is certainly somewhat complicated because of the non-local character of the operator within 
brackets. 

In the present context, however, there is a much simpler alternative which is equally well 
motivated. Rather than keeping a constant source and looking for the field it creates at a 
given point, we may ask the converse: What source J{x) will give a constant field (f){t) — >■ 1 



at late times? Classically, if the expansion rate is a constant, H = Hq, then a constant 
source will produce an asymptotically constant field J/{2Hq) at t ^ oo. In the 

quasi-de Sitter limit (p ^ 1), a constant field — >■ 1 is caused by a source of the form 
J — > 2H^{t). In this sense, the source which is needed in order to keep — >■ 1 can be used 
as a measure of the local expansion rate. Needless to say, this measure reduces trivially to 
the curvature scalar 

2H^{t) ^ Jit) = (D + ^n^i = ^n. (31) 

This illustrates the fact that in a quasi-de Sitter phase we may adopt ^n/l2 as a local 
definition of the expansion rate, because the second term in (1271) is negligible. The expansion 
rate can change by a large amount in the course of time, but as long as it does so adiabatically, 
the curvature scalar will be a good tracer of H{t). 

Now, in the quantum theory, the metric fiuctuates. If we adopt the expectation value of 
the classical expression as our definition of H{t), we have 

H\t) ^^m. (32) 

Naively, in pure gravity with a cosmological constant, we may expect to have a relation of 
the form 

(7^(x)) = 4A , (33) 

which would readily imply the constancy of H{t), with no room for a secular screening. 
Intuitively, Eq. flHHl) seems to follow from the Heisenberg equation of motion, which the 
field operator is supposed to satisfy identically . Nevertheless, the definition of a physically 
meaningful (TZ), and the proof of an equation of the form (1331) . involves a number of subtleties 
related to gauge invariance and renormalization (see also [l3|, [M]). A full discussion of this 
point is postponed to the next Section. 

Before closing, one comment is in order about the classical back-reaction. Classically, in 
a fiat FRW universe, we have 

;^ = -^(l + ^efr), (34) 

where Weff = p/p is the ratio of pressure p to energy density p. As mentioned at the 
beginning of this section, the classical back-reaction due to short-wavelength gravitational 
waves modifies the expansion law like a usual radiation field, with Wrad = 1/3. If the density 
in radiation is comparable to the cosmological term, then (1 + Wcs) will not be small, and 
from (!7r|) . 71 will not be a good tracer of H(t). At the classical level, 71 is constant, 
and hence it is completely insensitive to the classical back-reaction effect. In this respect, 
the observable originally proposed in Ref. given in Eq. ( IHUj) . cannot sense the traceless 
component of the energy momentum tensor, either. Clearly, (p =const. is one of the solutions 
when we assume J =const. Although there are various solutions for (p for a constant J, this 
variety is due to the degrees of freedom of the initial conditions, which are supposed to be 
irrelevant at late times. The study of alternative observables which are, at least, sensitive 
to the back-reaction effect caused by classical gravitational waves, is postponed for future 
work. 



IV. IS THERE A SECULAR SCREENING? 



As discussed above, the Ricci scalar 7?. is a good indicator of the adiabatic evolution of 
the expansion rate in a quasi-de Sitter phase. Here, we will discuss the calculation of its 
renormalized expectation value in the theory of pure gravity. For this purpose, it will be 
quite important to work with quantities which are invariant under diffeomorphisms. When 
we consider the gauge transformation 

^ x'^ = - , (35) 

TZ{x) transforms like TZ{x) — TZ{x) ~ TZ{x) + TZ^^x^- this sense TZ is not invariant. 
However, we do not really need to measure the value of curvature at a specified point 
in spacetime. Rather, it will be sufficient for our purposes to consider its value smeared 
over a certain sample volume. Let us introduce a window function W{x) which by definition 
transforms as a spacetime scalar. Then, for any spacetime scalar operator 0{x), the integral 
/ d'^x W{x)0{x) is manifestly gauge invariant.^ It will be useful to introduce the 
following notation for the expectation value of this quantity: 



{O}^ = ( phys 



j d^x^/^ Wix)0{x) 



phys 



A precise definition of the arbitrary physical state |phys) will be given below. 
The basic goal of this Section is to show that the equation 

(7^ren(a;))ly = 4A (1 (36) 



ren 



holds for any scalar window function W. The definition of the operators 7^ren(a;) and 1 
appearing in fl36|) requires explanation. In the path integral approach, we can calculate the 
n-point functions {h^y{xi) ■ ■ ■ h^y{xn)) to arbitrary loop order. All divergences in this calcu- 
lation can be reabsorbed by diagrams involving the vertices generated by the counterterms 
in Set- Still, such renormalized n-point functions will not be free from divergences in the 
coincidence limit, when two or more of the n points are brought to sit on top of each other. 
Since ^J—gTZ contains the coincidence limit of n-point functions of h^^, the counterterms in 
Set will fail to render a finite expectation value for ^/^TZ. This situation, of course, is not 
specific to gravity, and the problem is remedied once we introduce a probe field which cou- 
ples to the composite operator of interest. This will allow us to define a suitable regularized 
operator ^/^T^ren whose renormalized expectation value is finite. 

It is instructive to start by considering the simpler example of an interacting scalar field 
ij) in Minkowski spacetime [Sj. In this case, the two point function {ip{x)ip{x')) is finite 
after renormalization, but its coincidence limit {ip'^{x)) is still divergent. In the context 
of a single free scalar field, there is no counterterm to renormalize the value of (ip'^lx)). 
On the other hand, we can only measure this seemingly divergent quantity through some 
interaction. Let us therefore introduce a coupling to a probe scalar field via the interaction 



6 Note that W{x) = W^{x) = S'-^\x'' - )/ ^—g{x), is not a suitable window function. Although this 
is a scalar with respect to the transformation of x, it is also so with respect to the transformation of z. 
The bi-scalar transforms as SWz{x) = x^{x){dWz{x) / dx'^) + x^{z){dWz{x) / dz^), and for that reason 
J d^x y/^TZ{x)Wz{x) = 7?.(z) is not gauge invariant (in agreement with the discussion below Eq. ([55]) ). 



Lagrangian — Now, we can "measure" X{tlj'^{x)) as a contribution to the mass of 
the probe field 0. Here, the probe field is treated as classical, meaning that we neglect 
all the loop diagrams containing its propagator. We also assume that its amplitude is 
infinitesimally small. Now, g{ijj'^{x)) can be renormalized because the divergence in {ip^lx)) 
can be absorbed by a mass counterterm of the 0-field. Hence, we have found a regularized 
operator Xip^^^{x) — Xip'^{x) + Srn^ whose renormalized expectation value 

is finite by virtue of the probe field counterterm 5jC^ — — (5m^)0^/2. 

The same argument works for y/—gTZ. We consider a probe massless scalar field with 
the curvature coupling as we discussed in the preceding section. The action we add is 

S^ + Stt = ~J d'x^ + + SMligW) , (37) 

where the mass counterterm t^A^^f^f] is made up of curvature invariants, and will be further 
specified below. We may now define 

e7^'^- =^n + 6Ml, (38) 

whose renormalized expectation value may be thought of as the local value of the mass of 
the (f) field 

".L-nlW - ^^Si^ • (39) 



Since the volume {l)w = (phys | J ci'^a; l/r(a;)| phys) contains polynomials of /i^j, 
through y/—g, this quantity is also divergent. Hence, we have to renormalize this expression 
by adding a counterterm 5yoi, i.e. (F'^'^)iv- = (1 + Svoi)w- To be more precise, in order to 
renormalize the volume, we need to add another probe field to measure it. For example, 
we can consider a scalar field with ^ = as a probe. In this case the renormalized volume 
integral of its mass will measure the renormalized volume. 

The value of m^^^j.^^^-, [W] will change depending on the choice of the finite part of countert- 
erms. However, if there is a choice of counterterms in which the relation ^tt^^^-j.^^^^ [VT] = 4^A is 
maintained independently of the window function W, it is such renormalization conditions 
that are natural and appropriate for the theory that we are considering. Finite renormal- 
ization of local counterterms will correspond to introducing new interactions between the 
probe field and gravity, different from the original curvature coupling term. We shall not 
pursue the consideration of such interactions here. They correspond to higher order irrele- 
vant operators (which arc not expected to lead to infrared effects of the sort we are interested 
in). Thus, the basic question is whether we can choose local counterterms which make the 
renormalized value of iri^^j.^^^ to be constant. 



The key identity is 



f f d^^^.^^w{x)g,,{x)e'' 

J 6=0 J drinAx) 



-(7^-4A)+ 9^^iS9f+FP + S^^) \ , (40) 



where we assumed that the window function W{x) is independent of hf^„ = g^^ — gj^J . gj^J 
is the background metric which can be different from the de Sitter one as long as it solves 
the Einstein equations. The first equality follows from functional integration by parts. 
In the second equality, we have dropped the term —16i6{0) which arises from functional 
differentiation 5gfj,u{x) /5h^y{x). It is clear that this naively divergent term is local, and can 
be grouped together with Set- In fact, such functional derivative vanishes in dimensional 
regularization. All the variables in the path integral are to be understood as (+)-fields and 
the integral over (— )-fields has been abbreviated, since the CTP formalism is not essential 
for the current discussion. 

We submit that the appropriate choice of 5A^^ which implements our renormalization 
condition is given by 

8Ml = 2<4^^ + 4eA5™i . (41) 



Here, is the counterterm action for the theory of pure gravity, without the field. Thus, 
5A1^ is local as long as Set is so. The other part 4A(5voi is also local. Then, substituting in 
Eq. (l38l) . we have 



{n''^^)w = 4A (l^-)v^ - 2k ( 9>^^^yf+^P \ , (42) 



-9 /VK 

where we have used Eq. ( HOl) . 

The remaining task is to show that the second term in the right hand side of (H2l) vanishes 
when the expectation value is taken for physical states. To show this, the essential point 
is to understand what is meant by the physical state. Since is set to 0, we neglect it 
completely in the following discussion. It will be very convenient for our purposes to follow 
the standard construction for gauge fixing based on the BRST invariance p^]. The gauge 

transformation changes h^u{x) h^i,{x) = h^y{x) + 6h^^{x) with 

Sh^^{x) = ^^pX% + 9upX'',^, + 9tiu,pX'' ■ (43) 

The BRST transformation 5b oi h^v is obtained by simply replacing with a Grassmanian 
field in Eq. ( H3l) . The BRST transformation of is determined by requiring the nilpotency 
of the BRST transformation, 5\h^i, = 0. Different from the usual gauge theory, this equation 
does not determine 6bc'^ locally. The obtained equation contains derivatives of 6bc^- Hence, 
we do not give an explicit expression for 5bc^, which is not required below. We add the 
anti-ghost field and its BRST transformation introduces i^'^-field as 6bc^ = iB^. From 
the requirement of nilpotency of the BRST transformation, we have 5bB^ = 0. After these 
preparations, for an arbitrary gauge fixing function F^lhap], the gauge fixing term and the 
Faddeev-Popov ghost term are simultaneously given by 



Sgf + FP — I d X Cgf + FP , 

with 

^gf+FP = —i^B 



1 



5^ (f, + l-aB,] - i{5BKp)^-d' , (44) 

V ^ / 0hal3 



a/3 

where indices in Cgf+pp are raised and lowered by using the background metric g^^J . Since 
F'^ilha/s] may contain differentiation of h^fs, 6F^/6ha/3 is understood as the derivative operator 
that is obtained by the usual variational principle. In the present case it acts on c'^. For 
simplicity, we assume that F^[ha/3] is linear in h^/s. Hence, 6F^/6ha/3 is an operator solely 
written in terms of the background quantities. 

Let us now consider the physical observables and the physical states. In the BRST 
formalism, observables are BRST invariant quantities. This corresponds to the usual notion 
of gauge invariant variables such as the Bardeen parameter at the linear order. We should 
note that 6bs{x) ^ for a scalar s{x), which is the reason why we had to introduce a 
window function W{x) to evaluate the expectation value of TZ. An observable O satisfies 
[Qb,0} = 0, where Qb is the BRST charge defined in such a way that 6b* = [Qb,*}- 
Correspondingly, physical states are also required to be BRST invariant. Hence, they must 
satisfy 

g^lphys) = . 

Therefore, for physical states, any operator that can be written in the "exact" form [Qb, *} 
has vanishing expectation value. 

Then, using 6bW{x) = c'^dnW{x), which is the standard transformation rule for any 
scalar quantity, it is straightforward to show that 

S f 

d'^xWg^^-^ / d'^x Cgf+FP 
Oh,,,, J 




[ d'xw(g„,^B"-t{5Bh 



[iv riper 

Qb,-i [ d'xWg^,,^c^ \ . (45) 
J bh^„ J 

Hence, the contribution from Sgj+pp vanishes when the expectation value is calculated for 
physical states. This finally establishes our claim that we can choose local counterterms 
such that {7lren{x))w = 4A (Iren)vK holds for an arbitrary scalar window function W{x). 
This simply means that TZrenix), as measured by its effect on a probe scalar field, stays a 
constant over the entire space-time. 



V. CONCLUSION 

A secular screening of the cosmological constant by infrared quantum effects would rep- 
resent a very spectacular phenomenon in low energy quantum gravity. In this note, we have 
reanalyzed the issue of gauge invariance in the definition of the expansion rate H{t) which 
was used in the original analysis of this problem ^ (see Eq. ([7])). 

We have shown that such definition is only invariant under c-number gauge transforma- 
tions, but not under generic changes of the gauge fixing term. Such changes correspond to 
gauge transformations where the gauge parameter x'^ depends on the operator h^,,. Because 
of that, they introduce arbitrary time dependence in the expansion rate H[t) as defined in 



Eq. ([7]). Hence, the interpretation of the results in Ref. [2|] as a physical screening of A seems 
very questionable. 

A truly gauge invariant definition of H{t) was introduced in Ref. This definition 
was motivated on physical grounds as follows. A constant source J in a quasi-de Sitter 
universe coupled to a conformal scalar field will produce a field The amplitude of a 
free conformal scalar in quasi-de Sitter decays with time like the inverse of the scale factor. 
Hence, the late time behaviour of (j){t) is dominated by the contribution of the source during 
the last e-folding time, and is therefore proportional to the surface of a sphere of Hubble 
size 

(f){t) oc JH-\t) . 

It was proposed in that such auxiliary field be used as a measure of the local expansion 
rate. The field (p is given by the inverse of the perturbed wave operator acting on the 
constant source. This is a non-local and rather cumbersome expression to deal with in the 
quantum theory. On the other hand, we have argued that there is an alternative definition 
which is equally useful if we wish to monitor an adiabatic change in the expansion rate 
(such as the one which would be suggested by Eq. ([H])). Indeed the curvature scalar 71 is 
proportional to H'^{t) plus corrections of order H which are negligible in the adiabatic limit. 
So the question is whether the value of this scalar (or a suitable smearing of it) can change 
in the course of time. Classically, for the system of pure gravity coupled to a cosmological 
constant, this is impossible. By using the path integral approach, we confirmed that this 
conclusion is not altered when we take into account the subtleties associated with gauge 
invariance and renormalization. Therefore, according to this definition, we find no evidence 
of a secular screening of the cosmological constant, to all orders in perturbation theory. 

It should be stressed that these arguments apply only to the case of pure gravity with a 
cosmological constant, and they do not exclude the possibility of interesting infrared effects 
in theories with a different field content 0, 0, @, lH or due to non-perturbative effects 12 



Our considerations focused on the renormalized expectation value of the scalar curvature 
TZ, which is insensitive to the classical back-reaction effect due to a bath of gravitons. In 
future work, we would like to examine different gauge invariant indicators of the expansion 
rate fis'], which give a non- vanishing result depending on the choice of the initial state. 

Note added: 

After this paper was submitted to the archives, Tsamis and Woodard wrote a reply to 
it 1^, expressing some points of view which we do not share. 



First, they claim that we did not show that the renormalized Ricci scalar is constant, 
and that our Eq. (IHUl) is completely consistent with screening. The observable we calculate 
is the expectation value of the integral of the Ricci scalar over a region of space-time. This 
operator is divergent, and so we define the corresponding renormalized operator by standard 
techniques. We show that this agrees with the expectation value of the integral of a constant, 
over the same region. The equality holds order by order in the loop expansion. The region 
of space-time is itself arbitrary, as long as the same one is used on both sides of the equation. 
In our view, this means that the renormalized Ricci scalar, as measured by its effect on a 
probe scalar field, stays constant, in as precise a sense as can be made. Notice that this is 
exactly the condition for J to be constant with constant in Eq. fl26p . 

The authors of Ref. [l6| object that we use an external scalar window function W{x) 
in our definition of the gauge invariant operator (The reason for that is explained in our 



footnote [H]). This scalar is not constructed from the metric, and hence, the integral of 
^W{x)7l does not correspond to any observable of the theory: it depends on the particular 
choice of W{x). However the statement that our equality (IHUj) holds for any W{x) is, of 
course, independent of this choice, and hence it is a physically meaningful statement. They 
also object that even if we show that {7l'^^^)w = 4A {V^^)w , both sides of the equation 
can evolve secularly in the same way. Even if that were the case, this would not imply any 
secular evolution of their ratio, which is the quantity of our interest (for constant 0, the 
ratio is proportional to J at the classical level). They also claim, at the beginning of Section 
3 that our renormalization scheme is "peculiar". We disagree with this appreciation. What 
we do is standard renormalization in low energy effective theory. We do make a particular 
choice for the finite parts of the local counterterms which need to be subtracted. This 
is explained in detail in the paragraphs between our Eqs. ( l36l) and ( HOl) . The important 
point is that there is a choice of counterterms for which there is no secular screening of 
the renormalized operator. If a change in the local counterterms happened to give rise to 
some additional effect, then this would be an effect due to local physics (or, conceivably, 
to the secular evolution of the added higher-order local counter terms, although this seems 
unlikely), but it would be unrelated to the infrared secular evolution of 71. 

The authors of 1^ also purport that if we are allowed arbitrary subtractions in order 
to construct the renormalized operator TZ^en, then we could absorb in its definition things 
like the one loop effective potential of a scalar field. If so, they argue, we would reach the 
conclusion that T^ren stays constant even in a theory like " new inflation" , where the potential 
is due to one loop corrections. Of course this would not be correct, and it has nothing to do 
with the method we are using in the present paper. Arbitrary subtractions are simply not 
allowed. At each order in the loop expansion, we only allow as counterterms a finite number 
of higher dimension local operators, suppressed by corresponding powers of the Planck mass 
Mpi. The number of counterterms will be larger if we work at a higher order, because this is 
unavoidable in non-renormalizable theories. But these higher order counterterms can never 
absorb the lower order loop-corrections since the power of Mpi is different. 

We would agree that there are other observables one can look at. Our claim is that we 
see no evidence for a secular infrared screening in the observable we have analyzed. We 
should add that this is a better defined observable than the spatially averaged Hubble rate 
used in j2|. The authors of [l6| claim in Section 2 of their reply that gauge dependent 
quantities can have some physical content. While this is debatable, their discussion does 
not seem to warrant the preference of a gauge-dependent result over the gauge-invariant one 
we presented in this paper. 
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